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ABSTRACT 

Extrasolar planets are a natural extension of the interacting binaries towards the companions 
with very small masses and similar tools might be used to study them. Unfortunately, the 
generally accepted treatment of the reflection effect in interacting binaries is not very suitable to 
study cold objects irradiated by hot objects or extrasolar planets. 

The aim of this paper is to develop a simple model of the reflection effect which could be 
easily incorporated into the present codes for modeling of interacting binaries so that they can 
be used to study above mentioned objects. 

Our simple model of the reflection effect takes into account the reflection (scattering) , heating 
and heat redistribution over the surface of the irradiated object. The shape of the objects is 
described by the Roche potential and limb and gravity darkening can be taken into account. The 
orbital revolution and rotation of the planet with proper Doppler shifts for the scattered and 
thermal radiation are also accounted for. Subsequently, light-curves and/or spectra of exoplanets 
were modeled and the effects of the heat redistribution, limb darkening/brightening, (non-)grey 
albedo, and non-spherical shape were studied. Recent observations of HD189733b, WASP12b, 
and Wasp-19b were reproduced reasonably well. HD189733b has low Bond albedo and intense 
heat redistribution. Wasp-19b has low Bond albedo and low heat redistribution. 

We also calculate the exact Roche shapes and temperature distribution over the surface of all 
78 transiting extrasolar planets known so far. It is found that the departures from the sphere vary 
considerably within the sample. Departures of about 1% are common. In some cases (WASP- 
12b, WASP-19b, WASP-33b) departures can reach about 14, 12, and 8%, respectively. The mean 
temperatures of these planets also vary considerably from 300 K to 2600 K. The extreme cases 
are WASP-33b, WASP-12b, and WASP-18b with mean temperatures of about 2600, 2430, and 
2330 K, respectively. 

Subject headings: (stars:) binaries (including multiple): close — (stars:) binaries: eclipsing — stars: 
low-mass, brown dwarfs — (stars:) planetary systems 



Introduction 



There are sophisticated computer codes for 
calculating and inverting light curves or spec- 
tra of binary stars with various sh apes or ge- 
ometry including the Roche model (iLucvl 1 19681: 

Wilson fc Devinnevlll97lllWoodll97ltlMochnacki fc Dous 



Orosz fc Hauschildl 120001 IB radstreet fc Steelmar 
2002[|Pribulld2004lPavlovski et al.l2006tlTamuz"et all 
20061 ). The Wilson & Devinney (WD) code is most 
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1994; 



19721 : iRucinskil Il973t iHilll Il979t IPonner fc Etzell 



often used and is continuously being improved 
or m odified (|Kallrath et al.lll998l IPrsa fc Zwitter 
2005). The reflection effect is taken into account 
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4l$lrnost of these codes. 

The standard description of this effect is given 
in Wilson (Il990h . This effect is understood in the 
following way: A surface element of star A is irra- 
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diated by many surface elements of star B. A frac- 
tion of impingin g energy called bolometric albedo 
(lR.ucinskilll969h is converted into the heat which 
rises the local temperature and re-radiates the en- 
ergy on the day side of the object. Rest of the 
impinging energy is plunged into the star. Typ- 
ically, bolometric albedo is set to 1 for radiative 
and 0.5 for convective envelopes. An increase of 
the temperature on the day side of one object trig- 
gers a secondary reflection effect on the second ob- 
ject and one, two or several iterations are allowed 
to converge to a final state. Roche model, limb 
darkening and gravity darkening are taken into 
account. This model does a good job for many 
interacting binaries. 

However, I argue that the above mentioned 
model of the reflection effect should be revis- 
ited. There has been a lot of progress in the 
field since the time this reflection effect was de- 
veloped. New types of very cool objects such as 
brown dwarfs and extrasolar planets have been 
discovered. These new areas evolve rapidly and 
produce interesting results. In many respects, an 
extrasolar planet can be understood as an inter- 
acting binary with a small mass companion. It is 
thus an attractive idea to apply these new results 
from the extrasolar planets to interacting binaries 
and vice versa. For example, models of extrasolar 
planets can provide day-night heat redistribution 
and reflected light while models of interacting bi- 
naries can provide sophisticated Roche geometry. 

The standard model of the reflection effect faces 
several problems which prevent its application to 
cool objects irradiated by hot objects and to extra- 
solar planets. When it comes to very cold strongly 
irradiated objects a considerable amount of energy 
might be reflected off the surface. The above men- 
tioned reflection effect in the interacting binaries 
neglects this reflected light which can be essential, 
especially at the shorter wavelength. Well known 
example of such effect are the planets and moons 
in our Solar system in the visible light. This re- 
flected light bears the spectroscopic signatures of 
the hot irradiating star and is not converted into 
the heat and re-radiated. It is commonly taken 
into acc ount in models of hot- Jupiters or planets i n 
general ( Seager et al. 20001 Sudarskv et al.ll2005T) . 
The definition of the albedo in the planetary sci- 
ences is almost the opposite of its meaning in the 
interacting binaries. In planets, the albedo is a 



fraction of the impinging energy which is reflected 
off the object and is not absorbed and converted 
into the heat. Reflected and re-radiated photons 
may also have completely different Doppler shifts 
depending on the mutual velocities of the two ob- 
jects and observer. 

Moreover, some portion of the energy absorbed 
on the day side can be transferred to and irradi- 
ated from the night side. Various authors devel- 
oped different parametrization of thi s effect in con- 
necti o n with the extrasola r planets ( Guillot et al.l 
119961 IBurrows et ail 120061 [Cowan fc Agoll l2010l ). 
Calculat ions of atmosphere models of extrasolar 
planets (Irlubenv et al.1 12003 : Barman et al. 20051 
IBurrows et al.l 120081: iFortnev et al.l 120081 ) demon- 
strate that there is a deep temperature plateau 
which turns off the convection in the atmospheres. 
Convection and vertical energy transport operates 
only at very deep layers. Observations of hot- 
Jupiters indicate that only a very small fraction 
of the impinging radiation ( 10 ) could have bee n 
plunged into the object (IBurrows et al. 20071). 
Hydr o dynamical simulations dDobbs -Dixon fc Linl 
| 2008[ IShowman et al.l 120091 iMenou fc Rauscher 
2009) of the atmospheres of extrasolar planets 
indicate that there are very strong horizontal cur- 
rents and jets which can effectively redistribute 
and circulate the energy between the day and 
night side of the planet especially along the lines 
of constant latitudes. How much energy gets re- 
distributed to the night side depends mainly on 
the complex structure and dynamics of the surface 
layers. 

The above mentioned effects and findings must 
be taken into account in the new or revisited model 
of the reflection effect when modeling very cold 
components of interacting binaries and extrasolar 
planets. At the same time, many transiting exo- 
planets are so close to their host stars that their 
shape may depart from the sphere and be best 
described by the Roche model. 

The main purpose of this paper is to develop a 
new simple model of the reflection effect which 
would consider the reflected (scattered) light, 
heating (absorption of the light), and heat redistri- 
bution over the surface and which could be used to 
model the cold components of interacting binaries 
and extrasolar planets. The model was included 
into our code shell spec (jBudai fc Richardsl2004l: 
Budai et al.1 120051 ). This program was origi- 
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nally designed to calculate light-curves, spec- 
tra and images of interacting binaries immersed 
in a moving circumstellar environment which is 
optically thin. The code is freely available at 
|http : //www . ta3 . sk/^ budaj / shellspec . htmi| 
and might be used to study various effects ob- 
served or expected in the interacting binaries or 
extrasolar planets. A Fortran90 version of the 
code whic h also solves the invers e problem was 
created bv lTkachenko et ab ( 20ld ). 



2. Roche model 

In the shellspec code the Roche model serves 
as a boundary condition for the radiative transfer 
in the circumstellar matter. Both objects, star 
and companion, may have shapes according to 
the Roche model for detached or contact systems. 
D escriptions of the Roche model can be found 
in lKopall (|l959t ). IPlavec fc Kratochviil (|l964l ). and 
many other papers and books. Let us assume a 
Cartesian coordinate system (x,y,z) centered on 
one of the stars (labeled as 1) such that the com- 
panion (labeled as 2) is at (1,0,0) and revolves 
around the z axis in the direction of positive y 
axis. Let the mass ratio, q, always be m^/m-i or 
'companion/star' and q < 1 will indicate the com- 
panion is lighter while q > 1 means the central star 
is lighter. Then, the normalized Roche potential, 
C, is expressed as: 



2q 



(l + q)n (l + q)r 2 



1 



C(x,y,z) = 

2 

1 \ , 2 

+ y 



where 



y 



(i) 



(2) 



y/{x - l) 2 + y 2 + z' 2 . 

The Roche surface of a detached component is de- 
fined as an equipotential surface C s = C(x s , y Sl z s ) 
passing through the sub-stellar point (x s ,y s ,z s ) 
(point on the surface of the star in between the 
stars, < x s < l,y a = z s — 0) which is localized 
by the 'fill-in' parameter /* < 1. We define this 
by: 



f i — %s i ^Ixt 



h = (1 - 0/(1 - L lx ) (3) 



for the primary and the secondary, respectively. 
L\ x is the x coordinate of the LI point Ll(Li x , 0, 0). 



Derivatives of the Roche potential are: 
dC 



dC 
dx 



Cy — -p— , C z 

dy 



dC_ 

dz 



(4) 



Gravity darkening of the non spherical objects 
is taken into account by varying the surface tem- 
perature according to the von Zeipel's law: 



T/T p 



(9/9p) 



(5) 



where g is the normalized surface gravity, /3 is the 
gravity darkening coefficient, T p ,g p are the tem- 
perature and gravity at the rotation pole. The 
normalized gravity vector is g = (C Xl C y , C z ) and: 



g= JC 2 +C 2 +Cl 



(6) 



Notice, that there is an imminent singularity in the 
calculations in the vicinity of LI, L2 points since 
gravity falls to zero which drags temperatures (a 
denominator in many equations) to zero. We avoid 
the problem by setting the lowest possible value of 

g/g P = io- 4 . 

Limb darkening is taken into account using the 
quadratic limb darkening law: 

1(0) = 1(0) fLD (7) 

Ild = 1 -ui(l - cos 6) - u 2 (l -cos df (8) 

and by calculating the cosine of the angle 9 be- 
tween the line of sight unit vector n = (n x , n y ,n z ) 
and a normal to the surface: 

„ , n x C x + n y C y + n z C z 
cos9 = -n.g/g = y I — . (9) 



3. Irradiation and heat redistribution 

In this section we will describe our treatment 
of the mutual irradiation of the objects. It can be 
applied to both objects but we will neglect multi- 
ple reflections between the two objects since this 
is not essential if one of them is much less lumi- 
nous. We will distinguish between three separate 
processes: reflection of the light off the object 
(or scattering which does not produce any heating 
of the irradiated surface); heating of the irradi- 
ated surface (day side) by the absorbed light; and 
subsequent heat redistribution over the entire 
surface of the object. Let's assume that the day 
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side of a planet is irradiated by the star then the 
impinging flux at a location r from star at r* is: 



F 

7 7 



cos S- 



(r- 



(10) 



where i?*,T* are radius and effective temperature 
of the star, 8 is an irradiating angle which is the 
zenith distance of the center of the star as seen 
from the surface of the planet: 



cos S 



(r - r*).g 



|r-r*|p 



(ii) 



Now, let's define two local parameters Ab(oc,/3), 
P r (a,j3) on the day side of the planet (irradiated 
side of an object) where a,/3 are the longitude 
and latitude, respectively, measured from the sub- 
stellar point. P r will be the local heat redistri- 
bution parameter and Ab will be the local Bond 
albedo of the surface. Consequently AbF^ is flux 
immediately reflected off the surface, (1 — As)Fi r 
will be a fraction of the irradiating energy which 
is converted into the heat, P r (l — Ag)Fi r will be 
a part of the latter which is redistributed over the 
day and night side of the object, while the remain- 
ing part, (1 — P r )(l — A B )F ir , will heat the local 
area. Then the energy conservation for the day- 
night heat transport can be written as: 



P r (l - A B )F ir dS day = // aTl n dS day+ 



\-night 

(12) 

Let's assume, of simplicity (or in the absence of 
a better approximation), that Ab(cx, 0),P r (a, /3) 
are constant and the heat is homogeneously redis- 
tributed over the day and night sides so that the 
surface has a constant temperature T . Then, 

P r (l — Ab) JJ F ir dSday/ JJ dSday+night = 0"Tq . 

(13) 



1 Note that our local P r should not be conf used with the 
global P n parameter of Burrows ct al. (2006). P n is a frac- 
tion of the impinging stellar radiation which is transferred 
to and reradiated from the night side, P n «i P r (l — Ab)/2. 
P n is from the interval 0-0.5 while P r runs from to 1. P r 
is a direct indicator of the heat redistribution while the P n 
parameter only reflects a combination of the heat redistri- 
bution and Bond albedo. 



In case the planet has a spherical shape and is far 
from the star this reduces to: 



r n 4 = 



i 



A B )§Tl 



(14) 



Let's explore another case, namely that the hori- 
zontal circulation on the planet along the lines of 
constant latitude is so strong that it dominates the 
day-night heat transport and the equilibrium sur- 
face temperature, 7i, will be a function of latitude. 
In this case and the assumptions above (spherical 
planet far from the star), one can consider an en- 
ergy conservation for a fixed latitude: 

P r (l-A B ) J F„ cos((3)da = <jT^{(3) J cos((3)da 

(15) 

where 

cos(<5) = cos(a) cos(/3). (16) 

The solution is that the temperature depends on 
the fourth root of cos(/3): 

r x 4 (/3) = i Pr (l-A B )^cos(P) = 1t 4 cos(/3). 

(17) 

The study of these two extreme cases lead us to 
suggest a heat redistribution model in which the 
day-night heat transport is a linear combination of 
the two cases mentioned above. Namely, we will 
express the surface temperature in the following 
way: 

TL(P)= T o[Pa + P b cosm (18) 

where P a ,Pb are the 'zonal temperature redistri- 
bution parameters'. P a =< 0, 1 > and Pb is to be 
determined from Eq[l2] so that the total energy 
budget is conserved. From Eq. 1121131 and US] we 
obtain: 



Pb = (1 - Pa) 



J dS da 



y-\-night 



J COS ft 'dS 'day +night 



(19) 



It can be shown that in case of a spherical planet 
far from the star 



P b =-(l-Pa) 
7T 



(20) 



and Pb —< 0, 4/tt >. Notice, that P a is a measure 
of the effectiveness of the homogeneous tempera- 
ture distribution over the surface versus the zonal 
distribution. It is intimately linked with the ef- 
fectiveness of the heat flows along the meridians 
versus parallels. 
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Finally, the temperature distribution on the 
surface of the irradiated planet will be: 



rp4 rp4 



T, 



dn 



old 



(21) 



where Tf T = (1 - P r )(l - A B )F„/a on the day 
side, Tf r — on the night side, and T Q id is the 
prior temperature distribution over the surface in 
the absence of the irradiation (including the grav- 
ity darkening etc.). It should be noted that im- 
posing the external irradiation on one side of the 
object c an alter the origina l temperature distri- 
bution. iBudai et al. (2009) argue that the core 
cooling rates (T id) from the day and night side of 
a strongly irradiated planet may not be the same 
and that the difference depends on several impor- 
tant parameters, such as the effectiveness and the 
depth where the day-night heat transport occurs, 
the stellar irradiation flux, and vertical redistri- 
bution of the opacities, atmospheric abundances, 
and/or presence of the stratospheres. 

Once we know the temperature distribution 
over the surface, one can approximate the mono- 
chromatic flux from the surface as being composed 
of two parts: 



pre fleet 



F — F 



F 



'tire 



(22) 



Reflection depends on the surface albedo A v and 
has to take into account the mutual velocities: 



preflect _ a jp . 
± ii **-v 2 ± v\ ir 



(23) 

i? 2 

F Vlir = cos 5- ^F* (24) 

(r — r*) z 

where F* ± is properly shifted flux emerging from 
the surface of the irradiating star. The Doppler 
shifts are the following 



y-2 



V2 = 



(r-r*).(v- v*) 



(25) 



(26) 



(27) 



where v is the velocity field vector at the given 
point on the irradiated surface specified by the 
vector r, v* is the velocity of the center of mass 
of the irradiating star, and z coordinate points to 



the observer. @ To calculate the reflected intensity 
we assume that the reflection is isotropic in which 
case: 

jreflect = preflect ^ ( 2 8) 

Finally, F* hermal can be approximated by the flux 
emerging from the non-irradiated model atmo- 
sphere with the effective temperature equal to the 
surface temperature of the irradiated object given 
by EqUTJ 



= F„ 2 {T e}s =T) 



and the associated intensity is given by: 



thermal 



4(0) 



thermal 



LD 



^thermal 



(29) 



(30) 



(31) 



7r(l - Ml/3 - u 2 /6) 

In this way we fully include into account the mu- 
tual velocities of the objects and observer, the ro- 
tation of the reflecting object, its limb darkening 
but neglect the rotation of the irradiating object. 
In some cases the later can be easily taken into ac- 
count by feeding the code with the precalculated 
rotationally broadened spectrum F*. 

Local Bond albedo used here is a weighted 
monochromatic albedo A v : 



A B = 



\A v Fldv 
jF*dv • 



(32) 



One has to keep in mind that our albedo refers 
to the reflected light only (not to the absorbed 
and re-radiated light). In case the irradiated ob- 
ject is very cold compared to the irradiating ob- 
ject there is a clear distinction between its ther- 
mal radiation and reflected radiation. However, 
if the two objects have comparable temperatures 
it is almost impossible to distinguish whether a 
particular photon was scattered or absorbed and 
re-radiated. In this case it is still possible to use 
our formalism and e.g. approximate the albedo 
by the single scattering albedo to cope with the 
problem. 



2 Note that generally Doppler shifts in the scattered and 
thermal radiation are not the same and may not be trivial. 
More detailed treatement would require high resolution ra- 
diative transfer in the irradiated atmospheres for a set of 
radial planet-star velocities. 
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It might be convenient to define the mean tem- 
perature of the whole distorted object: 



T 



_ J T 4 dS 



(33) 



Note that it docs not takes into account the re- 
flected light and should not be confused with the 
effective temperature or brightness temperature. 

4. Application to extrasolar planets 

Extrasolar planets are an attractive test-bed 
for such models. They are a natural exten- 
sion of the interacting binaries towards cool and 
small mass companions. Their radiation is gov- 
erned by the stellar insolation. Very sophisti- 
cated atmosph ere models of exoplanets are be- 



catea atmospn ere moaeis a t exop 
ing developed rtSeager et al. 2000t 



2003 



Fortney et al 



Barman et al. 
20081) 



20051 : 



Hubenv et al 



Burrows et al.1 



v et ai. 
.1 120081: 



At the same time, mod- 
els with various simplifications are often ver 
various 



en very 

useful to study variou s effects (jHansenl 2008; 
Kipping fc Tinettil [200l iMadhusudhan fc Seager 
20091 : ICowan fc Agolll2010l ). 



4.1. Shapes of the transiting extrasolar 
planets 

Our knowledge about the size of exoplanets 
comes mainly from the observations of transits. 
A majority exoplanet studies assume that the 
exoplanets have a spherical shape. The preci- 
sion of the photometric measurements is advanc- 
ing revealing more and m ore details. Recently, 
Kipping fc Tinettil (|2009l ) attempted to model 
the light curve and transit assuming the spher- 
ical shape but taking into account the intrin- 
sic night side emission of the planet during the 
transit. They found that in the i nfrared region 



the eff e ct is of the order o f 10 . Seager fc Hui 



(2002), iBarnes fc Fortnevl (2003) C arter fc Winn 



201C ) assumed the sha pe of the rotational ellip- 



soid. Welsh et al. I (l2010h used the Roche potential 
approximation to describe the shape of the par- 
ent star and to explain the ellipsoidal variations 
of HAT-P-7. Nevertheless, many transiting exo- 
planets are so close to their host stars that their 
shape is best described by the Roche model. Very 
recently, iLi et all (|2010h developed a very sophis- 
ticated model for the tidally distorted exoplanet 
WASP-12b. 



Our code was applied to transiting extrasolar 
planets and exact Roche model shape of all the 
transiting exoplanets was calculated. Circular or- 
bit with the radius equal to the semi-major axis 
and synchronized rotation of the planet was as- 
sumed. The results are displayed in the Table 1. 
It lists R sub , R b ack, Rpoie, Rside, Ref f which are 
the radii at the sub-stellar point, anti-stellar point, 
rotation pole, on the side, and the effective radius 
of the planet, i.e. the radius of the sphere with the 
same volume as the Roche surface, respectively. 
Departures from the sphere are measured by the 
value of the R su b/ Rpoie- It was also assumed that 
the observed radius of the planet is equal to R S ide 
which has almost no effect on the relative pro- 
portions of most of the planets but may slightly 
underestimate the departures from the sphere for 
a few highly distorted planets. 

One can see that the departures from the spher- 
ical shape vary by several orders of magnitude. 
Departures of about 1% are common. About 
8% of planets (all of them have semi-major axis 
smaller than 0.03 AU) have departures larger than 
3%. These include: OGLE-TR-56b, WASP-4b, 
and CoRoT-lb. The extreme cases are WASP- 
12b, WASP-19b, and WASP-33b when the depar- 
tures can exceed 14%, 12% and 8%, respectively. 
This is clearly comparable to the precision of the 
planet radius determination or the transit radius 
effect. Moreover, this shape distortion would be 
even higher if one was to consider the eccentricity 
of the orbit and assume the periastron distance 
instead of the semi-major axis. The highest de- 
partures seems to be along the line joining the 
objects (the R su b/ Rpoie parameter) which would 
affect mainly the overall light-curve. During the 
transit event, it is mainly R S ide/ Rpoie parameter 
which is most relevant and this does not acquire 
such high values (about 2% in case of WASP-12b, 
WASP-19b, and WASP-33b). 

The observed radius of the planet should be as- 
sociated with the cross-section of the planet during 
the transit which is represented by R S id e , Rpoie- 
On the other hand, theoretical radius, e.g. from 
the evolutionary calculations, might be associated 
rather with the effective radius of the planet. One 
has to keep the findings above in mind when inter- 
preting the planet radius measurements or calcu- 
lation s ( Guillot fc Showman 20021: Burrows et al 



20071 : iFortnev et al.l 120071 LBaraffe et all 12008: 
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Leconte et al. 2010h . 

Apart from the shape of these planets, I cal- 
culated also a temperature distribution over the 
surface for all these exoplanets. This is rather 
cumbersome to tabulate and the result depends 
heavily on the free parameters of the model that 
is why I provide the reader with the code which 
can do that. In the Table 1 only the mean 
temperature of each transiting planet is listed. 
This was calculated assuming its Roche shape and 
A B = 0.1, P r = 0.5, P a = 0.5, T u = WOK. 
These mean temperatures vary considerably from 
300 to 2600K. The hottest planets are WASP-33b, 
WASP-12b, and WASP-18b with the mean tem- 
peratures of about 2600, 2430, and 2330 K, re- 
spectively. The separation of these extremely hot 
planets from the host star is smaller than 0.03 AU 
and temperatures of planets at a wider separation 
decrease accordingly. 

4.2. Light curves of the extrasolar planets 

In this section, we apply our model to the 
light-curves of extrasolar planets HD189733b and 
Wasp-12b. We will start with HD18 9733b. The 
planet pr operties were d e termin ed bv lWinn et al 



can proceed to a comparis o n be 



Now, we 

tween the observations of iKnutson et al 



(|2007h 



(|2007l ). IKnutson et all (|20071 ) obtained superb 



light curve of the planet which covers more than 
half of the orbit. It was observed at 8 micron in 
the infrared region. 

First, let us introduce the qualitative distri- 
bution of the temperature over the surface of 
HD189733b produced by our model. FigH] illus- 
trates the behaviour of the temperature for a rel- 
atively intense heat redistribution factor P r = 0.6 
and quite small zonal temperature redistribution 
parameter P a = 0.1 (i.e. intense flows along par- 
allels but not very intense heat flows along merid- 
ians). One can observe that this model produces 
hotter regions at the sub-stellar point and near 
the equator while cooler regions are at the poles. 
Fig|2] then illustrates the behaviour of the inten- 
sity. This is what would an observer see looking 
on the planet pole-on. Limb darkening was ap- 
plied to the model. One would see the dark polar 
regions, however, the hot sub-stellar point would 
not be the brightest since all hotter equatorial re- 
gions would be dim due to the limb darkening. 
Black body approximation was used to model the 
energy distribution. 



and models produced by SHELLSPEC shown in 
FigEl The synthetic light-curves are very sensitive 
to the heat redistribution parameter P r =< 0,1 >. 
Especially, at the phases close to the transit, 
when the night side of the planet is seen, the 
flux can change by the orders of magnitude. The 
best fit is obtained for P r = 0.5 — 0.8, which 
means that the surface heat redistribution is cru- 
cial and quite effective. T his finding is in agree - 
ment with the result s of IKnutson et all (12 007) . 
Burrows et al. ( 2008). iMadhusudhan fc Seagerl 



(|2009h . and lCowan fc Agoll (l201(lh . 

Albedo is an important parameter at all wave- 
lengths. The effect of the wavelength independent 
Bond albedo is illustrated in FigJH Our models 
indicate that Bond albedo of this planet is rela- 
tivel y small of about 0.1 which is in agreement 
with lCowan fc Agoll (|2010t ). Higher albedo reflects 
more light which could be observed especially at 
the shorter wavelengths. This reflected light might 
be relatively small in the IR region. However, 
higher albedo reduces the amount of energy ab- 
sorbed and redistributed over the surface. Surface 
temperatures are lower which manifests in lower 
fluxes at longer wavelengths. 

It is interesting to see that these kinds of light 
curves of transiting exoplanets are not very sen- 
sitive to the zonal temperature redistribution pa- 
rameter P a =< 0, 1 > associated with the effec- 
tiveness of the homogeneous heat transfer versus 
zonal transfer. It is illustrated in Fig|4]too. If the 
heat flows mainly along the equator and not along 
the meridians (P a = 0) and we view the planet 
almost edge-on then we observe a slight increase 
of flux at all phases, especially on the night side 
(compared to the homogeneous temperature dis- 
tribution, P a = 1). 

It is also interesting to observe a moderate ef- 
fect of the limb darkening on such light-curves of 
transiting exoplanets. Limb darkening is mani- 
fested mainly shortly before and after the transit 
and near the secondary eclipse. See Fig [5] with the 
comments and the description. It might be inter- 
esting to point out that extrasolar planets might 
show limb brightening at some wavelengths. This 
depends on the temperature gradient at a partic- 
ular depth probed by a certain wavelength. If the 
temperature is decreasing one would observe the 
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limb darkening. If it is increasing one would ob- 
serve limb brightening. There are two main ef- 
fects which could cause the temperature inversion. 
(1) Presence of species high in the atmosphere ef- 
fectively absorbing near the wavelength of max- 
imum of the stellar energy redistribution. This 
would cause the temperature inversion often called 
a stratosphere. (2) Heat redistribution between 
the day and the night side, especiall y if it occurs 
at deeper layers (|Burrows et alj|2008t ). This might 
significantly cool such layers on the day side pro- 
ducing a drop in the temperature and associated 
temperature inversion. This is the main motiva- 
tion why we study also the effect of limb bright- 
ening on the FigJS] and its effect is just the op- 
posite of the limb darkening. It might be worth 
mentioning that all the above mentioned calcu- 
lations took into account the proper inclination 
of the orbit. Depending how close to edge-on the 
planet's orbit is, the observed light curve may vary 
by 5% simply due to the observer's viewing angle 
(jCowan fc Ago]||2008h . 

However, there are cases when the zonal tem- 
perature redistribution parameter, P a , and limb 
darkening will be much more important. FigjB] 
illustrates the situation for a hypothetical planet 
at an inclination of 20 degrees. We used the pa- 
rameters for t he non-transi ti ng ex tr asolar planet 



HP 1 79949b (|Tinnev et al.1 l2001t iCowan et al 
120071 ) to feed the code with real numbers. How- 
ever, inclination, mass and radius are not known 
for non-transiting exoplanets. In this case the 
effects of the zonal temperature redistribution pa- 
rameter and limb darkening are indeed important 
and are comparable to the amplitude of the light 
curve. Another point worth making is that there 
will be a strong degeneracy between the radius, 
P r , and inclination since all of them affect the 
amplitude of the light-curve. Also there may be 
a degeneracy between P a , limb darkening and ra- 
dius since all of them affect the mean level of 
the light-curve. This justifies the complexity of 
our model with the heat redistribution and zonal 
temperature redistribution, and limb darkening 
parameters. Gravity darkening is also included 
but this will not be an issue for the extrasolar 
planets. If they are close to the star then their 
intrinsic radiation is negligible compared to the 
heavy external irradiation and if they are far from 
the star then their shape will be almost spherical 



(neglecting the rotation) and without the gravity 
darkening. 

Finally, let us study how this Roche shape 
might affect the light-curve of an extrasolar planet. 
This is illustrated FigJ7]on the example of WASP- 
12b at 8 micr on. The l ight v aries by about 10% as 
predicted by iLi et "all ( 2010| ). However, the most 
striking thing is the double humped light-curve in 
case of effective heat redistribution P r = 1. This 
shape of the light-curve gradually changes to a 
typical single humped shape for less intense heat 
redistribution (lower values of P r ). The transition 
to a single humped light-curve is color dependent 
and for shorter wavelengths occurs at higher P r . 
These calculations assumed the proper inclination 
of the orbit, zero albedo, zero limb darkening, and 
P a = 1. Assuming A B = 0, P r = 1, P a = 1, T old = 
100K and Roche shape we got the mean tempera- 
ture of the planet of 2470 K whi ch is in good agree- 
ment with 2516K obtained by iHebb et "ail (j2009l) 
who assumed the spherical shape. 

In the next step we calculated the light-curves 
of the planet at much shorter wavelength (0.9 mi- 
cron) and attempted to understand the secondary 
eclipse observations of Lopez-Morales et all J2009) . 
This is illustrated in Fig{8j Notice that at this 
shorter wavelength the light-curves cover a consid- 
erable range of values during the secondary eclipse. 
This is contrary to the longer wavelengths when 
light-curves have higher spred during the transit. 
There is a lot of degeneracy between the albedo 
and heat redistributi on in this case. For exam- 
ple the observation of lLopez-Morales et al. (2009) 
can be reproduced equally well by the model with 
Ab = 0,P r = 0.25 with the mean surface tem- 
perature 2470 K as well as by the model with 
Ab = 0.95, P r = 1 with the mean surface tem- 
perature of 1170K. Note that heat redistribution 
is not very important nor well constrained if the 
albedo is very high. It is also possible to fit the 
observations assuming that the surface is not re- 
flective and the temperature is homogeneous over 
the entire surface (Ab = 0,P r = l,P a — 1) and 
adopting T u — 2600-fT. This corresponds to the 
mean temperature of the planet of about 3000K 
(in this case the mean temperature is identical to 
the effective temperature and to the brightness 
temperature of the body). This is much higher 
than a day si de brightness temper a ture o f 2660 
K reported by lLopez-Morales et all (|2009h but it 
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is in per fect agreement with t he value of 2997 K 
given by ICowan fc Ago! (|2010l ) . 



4.3. Spectra of the extrasolar planets 

I would like to point out that this paper does 
not deal with the calculations of the local atmo- 
sphere models and intensity emerging from the lo- 
cal atmosphere models on the surface of the planet 
or star. However, in case grids of such mod- 
els and spectra are avail able (jHubenv fc Burrows 
2007 : Allard et al. 20031 ) they can be used as an 
input to our model and composite spectra of the 
star and/or planet can be calculated from optional 
view point. In this respect, our code is not re- 
stricted to the black body approximation. 

In this section we apply our model to WASP- 
19b and study its spectru m . Th is planet was 
discovered by iHebb et al.l (j2010h . The sec- 
ondary eclipse was de tected in the H-band by 
Anderson et al. I (l2010bh . The authors found the 
emission from the planet very strong and it 
presents an interesting puzzle for the models. This 
is a highly distorted exoplanet and we assumed 
a proper shape according to the Roche poten- 
tial, proper inclination, zero heat redistribution, 
homogeneous zonal temperature redistribution, 
zero limb and gravity darkening and calculated 
spectra during the secondary eclipse. It was also 
assumed that the star h as a spectrum given by 
Castelli fc Kurucd (12004 . As a first approxima- 
tion for the planet I used the black body approx- 
imation. The comparison of the data with our 
model is in FigJH Our black body model with 
zero albedo and zero heat redistribution P r = 
still under-predicts the radiation from the planet 
but it is within 2a from the measurements. 

This figure also illustrates the effect of the 
albedo on the spectrum. For example, rela- 
tively small frequency independent (grey) albedo 
of about 0.2 will strongly increase the planet radi- 
ation at the shorter wavelengths but decrease the 
fluxes from the planet at wavelengths longer than 
0.9 micron. This is because the non-zero albedo 
reduces the surface temperature of the planet but 
reflects the stellar light at shorter wavelengths. In 
the next step we assumed that the albedo is not 
grey but has the following color dependence 



4)(A/A r 



(1-^o) + ^o(A/A )t 



(34) 



where A$ is albedo at some reference wavelength 
Ao- We set 7 = — 4 to mimic the Rayleigh scat- 
tering with with A -4 dependence and adjust Aq 
to recover the same Bond albedo (0.2) as in the 
previous case. This case produces identical sur- 
face temperatures as in the grey albedo case but 
radiation below 0.6 micron is considerably higher 
than in the grey case, radiation above the 0.6 
micron is lower than in the grey case, and at 
still longer wavelength the relative difference be- 
tween the grey and color albedo calculations grad- 
ually decreases. Finally, I refrained from the 
black-body approximation for the planet and used 
the grid of non-irradiated atmosphere models of 
Hubenv fc Burrows! (|2007t) for T eff < 1800^ and 
Allard et al.1 (l2003h (BT-Settl) for T eff > 1800/4T. 
W ith this approach one c an fit the observations 
of Anderson et al.l (|2010bl ) reasonably well. The 
fit requires that the planet has low Bond albedo 
and low heat redistribution. The stratosphere, if 
present, should not have drastic impact on the 
spectrum since if the JHK band were reversed and 
were in absorption (instead of emission) it might 
be difficult to conserve the total flux and fit these 
observations simultaneously. 

5. Conclusions - Summary 

I argue that the reflection effect in the interact- 
ing binaries must be revisited in order to describe 
properly the radiation from the cool object irradi- 
ated by the hot object. Subsequently, a new model 
for the reflection effect was proposed and applied 
to an extreme case of an interacting binary - ex- 
trasolar planet. 

The new model introduces several free param- 
eters. Some of them are well known in the field of 
extrasolar planets. (1) Bond albedo Ab =< 0, 1 > 
which controls how much energy is reflected and 
how much is converted into the heat. Ab = 1 
means that all impinging energy is reflected off the 
surface and nothing gets converted into the heat. 
(2) Heat redistribution parameter P r =< 0,1 > 
which controls how much heat from a certain point 
is redistributed to other places and how much is re- 
radiated locally. P r — means that all absorbed 
heat is re-radiated locally and nothing gets trans- 
ported to other places. (3) Zonal temperature re- 
distribution parameter P a =< 0, 1 > is a mea- 
sure of the effectiveness of the homogeneous heat 
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redistribution over the surface versus the zonal 
distribution. P a = 1 means that the heat is ho- 
mogeneously redistributed over the surface while 
P a = means zonal redstribution and that the 
heat flows only along the parallels. Our model 
properly describes the shapes of the objects by 
means of the Roche potential and takes into ac- 
count gravity and limb darkening. At the same 
time it takes into account the orbital revolution, 
rotation and proper Doppler shifts in the scattered 
and thermal radiation. 

It is demonstrated on HD189733b that the 
light-curves of transiting extrasolar planets are 
mainly sensitive to the heat redistribution param- 
eter P r and not so sensitive to the zonal temper- 
ature redistribution parameter P a or to the limb 
darkening. However, light-curves of planets seen 
at very low inclinations are very sensitive also to 
the zonal temperature redistribution parameter P a 
as well as to the limb darkening. This planet has 
low Bond albedo and relatively intense heat redis- 
tribution. 

The effect of non-spherical shape on the light- 
curve can also be important. For highly distorted 
planets like WASP-12b this might cause a double 
humped curve with the amplitude of about 10% 
superposed on other types of variability. Observa- 
tions of this planet at 0.9 micron cannot constrain 
the Bond albedo or heat redistribution well. 

The effect of the grey and non-grey albedo on 
the spectrum of the Wasp-19b was studied as well. 
It was demonstrated that non-blackbody model 
with flu xes given by non-irradiate d atmosphere 
mode ls ( Hubenv fe Burrows 2007} Allard et al 



2003) can fit the observations of lAnderson et al 



( 2010bh reasonably well. The planet has low Bond 
albedo and low heat redistribution. 

We calculated the exact Roche sh apes of all 
curre ntly known transiting exoplanets (| Schneider 
1995) and found out that the departures from 
spherical symmetry may vary significantly. Depar- 
tures of the order of 1% are common, and can ex- 
ceed about 8% in most extreme cases like WASP- 
12b, WASP-19b, and WASP-33b. About 8% of 
transiting planets have departures more than 3% 
(all have semi-major axes smaller than 0.03 AU). 
Temperature redistribution over the surface of all 
these planets was also calculated. The mean 
temperatures of the planets also vary consider- 
ably from 300 K to 2600 K. The extreme cases 



are WASP-33b, WASP-12b, and WASP-18b with 
mean temperatures of about 2600, 2430, and 2330 
K, respectively. 

I thank Nick Cowan for valuable comments on 
the manuscript and Heather Knutson for sharing 
a copy of her data. The author acknowledges the 
support from the Marie Curie International Rein- 
tegration grant FP7-200297 and partial support 
from the VEGA grants 2/0078/10, 2/0074/09. 
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Table 1: Shapes of the transiting exoplanets. Columns are: a -semi- major axis in [AU], q -star/planet 
mass ratio, R su b -planet radius at the sub-stellar point (R su bflfl), Rback -planet radius at the anti-stellar 
point, Rpoie -planet radius at the rotation pole (0,0,i? po ; e ), Rside -planet radius at the side point (0,i? s id e ,0), 
(assumed equal to the planet radius determined from the transit), R e ff -effective radius of the planet, 
rr = Rsut/Rpoie -departure from the sphere, /, = R S ub/L\x -fill-in parameter of the Roche lobe, T -mean 
temperature of the planet in K. Radii are in units of Jupiter radius. Coordinates of the center of mass of 
the planet are (0,0,0). See the text for a more detailed information. 



GJ 1214b 

WASP-19b 

CoRoT-7b 

WASP-18b 

OGLE-TR-56b 

TrES-3 

WASP-12b 

OGLE-TR-113b 
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Fig. 1. — Example of a 2D projection map of the 
surface temperature of HD189733b calculated for 
P r = 0.6, P a = 0.1. x-axis points to the star, z-axis 
to the rotation pole. One can clearly see the hot 
region on the day side of the planet facing the star. 
The night side has a non zero temperature due to 
the efficient heat redistribution (P r — 0.6). The 
temperature at the polar region, (x,y) — (0,0), 
drops significantly since the heat is redistributed 
efficiently along parallels but much less effectively 
along the meridians (P a = 0.1). 
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Fig. 2. — Example of a 2D projection image of 
HD189733b at 8 micron seen pole-on. The sur- 
face intensity is in erg/s/Hz/cm 2 /sterad. It cor- 
responds to the temperate distribution from Fig|T] 
The limb darkening (m = 0.6, U2 = 0.2) was ap- 
plied to it. Notice that hot regions are now actu- 
ally dark due to the limb darkening. 





HD189733b . \ 






.' ? ¥y.*»*- . >» 


\ 




•y v^4v*.\v*.- / 

'•■ - : ; r^- \ 1 






\\ ■. • * s \\ // 1 Pr=0.8 

■ / / Pr=0.5 

\ \ \ ../ / / Pr=0.2 






\ \ '**..., / / Pr=0.1 

















Fig. 3.— Theoretical light-curves of HD189733b 
compa red with the observations of iKnutson et al. 
(|2007l ) at 8 micron (black points). Observe the 



strong dependence on the heat redistribution pa- 
rameter, P r , which could reach several orders of 
magnitude at phase zero. As = 0,P a = 1-0, and 
zero limb darkening were assumed. 




Fig. 4. — The effect of the Bond albedo and 
zonal temperature distribution on the light-curve 
of transiting exoplanet HD189733b at 8 micron. 
Higher Bond albedo reflects more light which is 
seen mainly at shorter wavelength. Consequently, 
less energy is available to be absorbed and redis- 
tributed over the surface. Temperatures are lower 
which means lower fluxes in the IR region. If 
P a = and the heat flows mainly along the par- 
allels then (compared to the homogeneous flows 
with P a = 1) one would detect slightly more 
light at all phases, especially on the night side 



Black points are the observations of lKnutson et al 
(|2007l ). P r — 0.6 and zero limb darkening were as- 
sumed. 
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Fig. 5. — The effect of the limb darken- 
ing/brightening on the theoretical light-curve of 
the transiting exoplanet HD189733b at 8 micron. 
Limb darkening has almost no effect at phase zero 
but much stronger effect shortly before and after 
the transit and near the secondary eclipse. This 
is because when we observe the night side with al- 
most constant temperature distribution the dark- 
ening at the limb is compensated by the bright- 
ening at the center to conserve the flux. How- 
ever, during the secondary eclipse the limb dark- 
ening suppresses the radiation from the cold re- 
gions and amplifies the radiation from the hot re- 
gions. The net effect is that planet is brighter 
at the phase 0.5. The opposite happens shortly 
before and after the transit. The possible limb 
brightening would have the oppo site behaviour. 



Black points are the observations of lKnutson et al 
(|2007t) . A B = 0.1, P r = 0.6, P a = 1.0 were as- 
sumed. 
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Fig. 6. — The effect of the zonal temperature re- 
distribution (P a parameter) and limb darkening on 
the light-curve of a hypothetical exoplanet seen at 
small inclination of 20 degrees. Contrary to the 
case with homogeneous temperature distribution 
(P a — 1), one would observe less flux from the 
models with the dominant east-west heat circula- 
tion (P a = 0) since the polar regions visible at 
this inclination are quite cool. Limb darkening 
would dim the hotter equatorial regions even fur- 
ther and reduce the observed flux, especially for 
models with the dominant east-west heat circula- 
tion (P a = 0). Ab = 0, P r = 0.5 were assumed. 
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Fig. 7. — Illustration of the effect of the non- 
spherical shape on the light-curve of WASP-12b 
at 8 microns. In case of effective heat redistribu- 
tion (P r = 1), one would observe a highly atypical 
double humped light-curve caused by the Roche 
shape. This resembles the well known ellipsoidal 
variations in interacting binaries. For lower val- 
ues of P r and less effective heat redistribution the 
light curve would acquire a more typical 'cosine' 
shape. Ab = 0, P a = 1 and zero limb darkening 
were assumed. 
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Fig. 8.— Light-curves o f WASP-12b at 0.9 
microns. Observations of lLopez-Morales et al 



(l2009h . jlack point, can be understood within sev- 
eral very different models. Low albedo and low 
heat redistribution model (As = 0,P r = 0.25) 
produces almost identical light-curve than high 
albedo models (Ab — 0.95, P r is not very impor- 
tant in this case) . Also a non- reflective model with 
homogeneous temperature distribution of about 
3000 K (A B = 0,P r = l,P a = 1) can reproduce 
the observations. Roche shape and zero limb dark- 
ening were assumed. 
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Fig. 9. — Planet to star flux ratio as a func- 
tion of wavelength for WASP-19b. Model (2) is 
a black-body model with zero albedo. Model (3) 
is a grey albedo model (Ab — 0.2). Model (4) 
is a non-grey albedo model (Rayleigh scattering) 
and the same Bond albedo, Ab — 0.2. Model 
(1) is a non black-bo dy model based on non- 
irradiated atmosph eres ( Hubenv fc Burrows1l2007 ; 
lAllard et al.l 2003|). Bla c k symb ol is a measure- 
ment of lAnderson et al.l (|2010bl ) in the H-band. 
Zero heat and homogeneous zonal temperature re- 
distribution (P r — 0,P a = 1) and zero limb dark- 
ening were assumed. 
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